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Abstract 


We construct and estimate a joint model of macroeconomic and yield curve 
dynamics. A small-scale rational expectations model describes the macroecon- 
omy. Bond yields are affine functions of the state variables of the macromodel, 
and are derived assuming absence of arbitrage opportunities and a flexible price 
of risk specification. While maintaining the tractability of the affine set-up, our 
approach provides a way to interpret yield dynamics in terms of macroeconomic 
fundamentals; time-varying risk premia, in particular, are associated with the 
fundamental sources of risk in the economy. In an application to German data, 
the model is able to capture the salient features of the term structure of interest 
rates and its forecasting performance is often superior to that of the best avail- 
able models based on latent factors. The model has also considerable success in 
accounting for features of the data that represent a puzzle for the expectations 
hypothesis. 


JEL classification: E43, E44, E47 
Keywords: Affine term-structure models, policy rules, new neo-classical 
synthesis 


Non-technical summary 


This paper aims at deepening our understanding of how macroeconomic fac- 
tors drive movements in the term structure of interest rates and how they affect the 
behavior of risk premia embedded in observed yields. While this has long been a 
topic on the agenda of both financial and macro economists, the focus of research 
as well as the methods used have often been different. On the one hand, financial 
economists have mainly focused on forecasting and pricing interest rate related se- 
curities. They have therefore developed powerful models based on the assumption 
of absence of arbitrage opportunities, but typically left unspecified the relationship 
between the term structure and other economic variables. Macro economists, on the 
other hand, have focused on understanding the relationship between interest rates, 
monetary policy and macroeconomic fundamentals. In so doing, however, they have 
typically not imposed restrictions to preclude arbitrage opportunities and they have 
often ignored the role of time-varying risk premia as an important component in 
explaining movements in yields over time. In other words, macro economists have 
typically relied on the “expectations hypothesis” (i.e. assumed zero or constant risk 
premia) in spite of its poor empirical record. 

This paper combines the two lines of research and presents a unified empiri- 
cal framework where a small structural model of the macro economy is combined 
with an arbitrage-free model of bond yields. The proposed model extends the term 
structure literature, since it shows how to derive bond prices using no-arbitrage 
conditions based on an explicit structural macroeconomic model, including both 
forward-looking and backward-looking elements. At the same time, we extend the 
macroeconomic literature by studying the term structure implications of a standard 
macro model within a dynamic no-arbitrage framework. The fact that we use a 
structural macroeconomic framework, rather than a reduced-form VAR. representa- 
tion of the data, is one of the main innovative features of our paper. One of the 
advantages of this approach is that it facilitates an economic interpretation of the 


results. 
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In an empirical application using German data, we show that there are syn- 
ergies to be exploited from modelling the term structure of interest rates using a 
combination of macro and finance approaches, and that this gives rise to sensible 
results. Notably, we show that our estimates of macroeconomic parameters, that are 
partly determined by the term-structure data, are consistent with those that would 
be estimated using only macroeconomic information. At the same time, our model’s 
explanatory power for the term-structure is comparable to that of state-of-the-art 
term-structure models based only on unobservable variables. 

We find that our model performs very well in terms of two evaluation criteria 
that we employ, namely the model’s out-of-sample yield forecasting performance and 
the ability of the model to account for deviations from the expectations hypothesis 
in the observed term structure. We argue that the model’s success is due to both the 
inclusion of macroeconomic variables in the information set and to the imposition of 
a large number of no-arbitrage and structural restrictions. The fact that we are able 
to account for deviations from the expectations hypothesis demonstrates that the 
dynamics of stochastic risk premia are important determinants of yield dynamics. 
Since risk premia in our model are driven by macroeconomic fundamentals, we con- 
clude that the dynamics of such premia can ultimately be attributed to underlying 


macroeconomic dynamics within a consistent framework. 
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1 Introduction 


Understanding the term structure of interest rates has long been a topic on the 
agenda of both financial and macro economists, albeit for different reasons. On the 
one hand, financial economists have mainly focused on forecasting and pricing inter- 
est rate related securities. They have therefore developed powerful models based on 
the assumption of absence of arbitrage opportunities, but typically left unspecified 
the relationship between the term structure and other economic variables. Macro 
economists, on the other hand, have focused on understanding the relationship be- 
tween interest rates, monetary policy and macroeconomic fundamentals. In so doing, 
however, they have typically relied on the “expectations hypothesis,” in spite of its 
poor empirical record. Combining these two lines of research seems fruitful, in that 
there are potential gains going both ways. 

This paper aims at presenting a unified empirical framework where a small 
structural model of the macro economy is combined with an arbitrage-free model 
of bond yields. We build on the work of Piazzesi (2003) and Ang and Piazzesi 
(2003), who introduce macroeconomic variables into the standard affine term struc- 
ture framework based on latent factors — e.g. Duffie and Kan (1996) and Dai and 
Singleton (2000). The main innovative feature of our paper is that we use a structural 
macroeconomic framework rather than starting from a reduced-form VAR represen- 
tation of the data. One of the advantages of this approach is to allow us to relax 
Ang and Piazzesi’s restriction that inflation and output be independent of the policy 
interest rate, thus facilitating an economic interpretation of the results. Our frame- 
work is similar in spirit to that in Wu (2002), who prices bonds within a calibrated 
rational expectations macro-model. The difference is that we estimate our model 
and allow a more empirically oriented specification of both the macro economy and 
the market price of risk.! 


Our estimation results, based on German data, show that macroeconomic fac- 


1A framework similar to ours is also employed in recent papers by Rudebusch and Wu (2003), 
who interpret latent term structure factors in terms of macroeconomic variables, and by Bekaert, 
Cho and Moreno (2003), who mix a structural macro framework with unobservable term structure 
factors. 
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tors affect the term-structure of interest rates in different ways. Monetary policy 
shocks have a marked impact on yields at short maturities, and a small effect at 
longer maturities. Inflation and output shocks mostly affect the curvature of the 
yield curve at medium-term maturities. Changes in the perceived inflation target 
have more lasting effects and tend to have a stronger impact on longer term yields. 

Our results also suggest that including macroeconomic variables in the infor- 
mation set helps to forecast yields. The out-of-sample forecasting performance of 
our model is superior to that of the best available affine term structure models for 
most maturities /horizons. 

Finally, we show that the risk premia generated by our model are sensible. First, 
the model can account for the features of the data which represent a puzzle for the 
expectations hypothesis, namely the finding of a negative and large — rather than 
positive and unit — coefficients obtained, for example by Campbell and Shiller (1991), 
in regressions of the yield change on the slope of the curve. Second, regressions based 
on risk-adjusted yields do, by and large, recover slope coefficients close to unity, i.e. 
the value consistent with the rational expectations hypothesis. 

The rest of the paper is organized as follows. Section 2 describes the main 
features of our general theoretical approach and then provides a brief overview of 
our estimation method. It also discusses the specific macroeconomic model which we 
employ in our empirical application. The estimation results, based on our application 
to German data is described in Section 3. Section 4 then discusses the forecasting 
performance of our model, compared to leading available alternatives. The ability of 


the model to solve the expectations puzzle is tested in Section 5. Section 6 concludes. 


2 The approach 


In recent years, the finance literature on the term structure of interest rates has made 
tremendous progress in a number of directions (see e.g. Dai and Singleton, 2003). 
Following the seminal paper by Duffie and Kan (1996), one of the most successful 
avenues of research has focused on models where the yields are affine functions of a 


vector of state variables. This literature, however, has typically not investigated the 
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connections between term structure and macroeconomic dynamics. In the rare cases 
in which macroeconomic variables—notably, the inflation rate have been included in 
estimated term-structure models, those variables have been modelled exogenously 
(e.g. Evans, 2003; Zaffaroni, 2001; Ang and Bekaert, 2004). The interactions be- 
tween macroeconomic and term structure dynamics have also been left unexplored 
in the macroeconomic literature, in spite of the fact that simple “policy rules” have 
often scored well in describing the dynamics of the short-term interest rate (e.g., 
Clarida, Galf and Gertler, 2000). 

An attempt to bridge this gap within an estimated, arbitrage-free framework 
has recently been made by Ang and Piazzesi (2003).? Those authors estimate a term 
structure model based on the assumption that the short term rate is affected partly 
by macroeconomic variables, as in the literature on simple monetary policy rules, 
and partly by unobservable factors, as in the affine term-structure literature. Ang 
and Piazzesi’s results suggest that macroeconomic variables have an important ex- 
planatory role for yields and that the inclusion of such variables in a term structure 
model can improve its one-step ahead forecasting performance. Nevertheless, unob- 
servable factors without a clear economic interpretation still play an important role 
in their model. Moreover, Ang and Piazzesi’s two-stage estimation method relies 
on the assumption that the short term interest rate does not affect macroeconomic 
variables. 

In order to redress these shortcomings, we construct a dynamic term structure 
model entirely based on macroeconomic factors, which allows for an explicit feedback 
from the short term (policy) rate to macroeconomic outcomes. The joint modelling 
of three key macroeconomic variables—namely, inflation, the output gap and the short 
term “policy” interest rate-should allow us to obtain a more accurate (endogenous) 
description of the dynamics of the short term rate. At the same time, our explicit 
modelling or risk premia should also help us in capturing the dynamics of the entire 
term-structure. 


In this section, we present our approach to model jointly the macroeconomy 


In related papers, Dewachter and Lyrio (2002) and Dewachter, Lyrio and Maes (2002) also 
estimate jointly a term structure model built on a continuous time VAR. 
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and the term structure. The main assumption we impose is that aggregate macro- 
economic relationships can be described using a linear framework. To motivate our 
approach, we start with an outline of the macroeconomic model that we use in our 
empirical analysis. We then cast this macro-model in a more general framework 
and show how to price bonds within such a framework based on the assumption of 


absence of arbitrage opportunities. 


2.1 A simple backward/forward looking macroeconomic model 


We rely on a structural macroeconomic model, whose choice is motivated by the fact 
that it could be derived from first principles. The model is certainly too stylized 
— for example in its ignoring foreign variables or the exchange rate — to provide 
a fully-satisfactory account of German macroeconomic dynamics. Nevertheless, it 
does include the minimal structure of a macroeconomic model proper. Our results 
in sections in Sections 4 and 5 suggest that such minimal structure does capture the 
central features of the dynamics of yields. 

The model of the economy includes just two equations which describe the evo- 


lution of inflation, m+, and the output gap, zr: 


Te = Hau [T1] + (dom Hr) ed gt + Et» 


Be = py yey + (1 — py) tea — Cp (re — Er [rey ]) + EF. 


The inflation equation implies that prices will be set as a markup on marginal 
cost, captured by the output gap term in the equation. The assumption of price 
stickiness generates the expected inflation term, while the lags capture inflation in- 
ertia. The output gap equation provides a description of the dynamics of aggregate 
demand, which is assumed to be affected by movements in the short term real inter- 
est rate. The forward looking term captures the intertemporal smoothing motives 
characterizing consumption, the main component of aggregate demand.? 

3 Both equations can be derived from first principles. More precisely, the inflation equation can be 
derived as the first order condition of the price-setting decision of firms acting in an environment with 


monopolistic competition. Monopolistic competition implies that prices will be set as a markup on 
marginal cost, which explains the presence of the output gap term in the equation. The assumption 
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The two equation above are often interpreted as appropriate to describe yearly 
data. Since we will employ monthly data in estimation, we recast the model at the 
monthly frequency using an approach similar to Rudebusch (2002). The equations 


that we will actually estimate are therefore 


12 3 
Tt = = 3 E [mii] + (1 — u) S but +ÔrTt + EF, (1) 
i=1 i=1 
M 12 3 
m = Lo 3 E; [ziy] + 1 he) X Gritti — Cr (re — Er [messi] +f, (2) 
1=1 i=1 


where all variables now are expressed at the monthly frequency, and where inflation 
is defined as m; = In P; —ln P;_12, where P, is the price level at t.4 The two equations 
include a forward-looking term capturing expectations over the next twelve months 
of inflation and output, respectively. The backward-looking components of the two 
equations are restricted to include only 3 lags of the dependent variable. This 
choice results in a more parsimonious empirical model. In the estimation, we impose 
Ur + (1 — pe) D3; Oni = 1, a version of the natural rate hypothesis. 


Finally, we need an assumption on how monetary policy is conducted in order 


of sticky prices generates the expected inflation term, as firms do not know when their prices 
will adjust next and therefore need do maximize the sum of current and expected future profits. 
The additional lagged inflation rate has been motivated through the assumption of partial price 
indexation (Christiano, Eichenbaum and Evans, 2001) or the presence of a set of firms that use a 
backward-looking rule of thumb to set prices (Galf and Gertler, 1999). The output gap equation 
can be derived from an intertemporal consumption Euler equation. The first term on the right-hand 
side is essentially Hall’s (1978) random walk hypothesis which states that consumption is equal to 
expected consumption tomorrow (in simple, closed-economy models, consumption equals output in 
equilibrium). This hypothesis is supplemented with two additional terms. First, a real interest rate 
(which Hall assumed to be constant) shifts the consumption profile such that a real rate increase 
tends to discourage current consumption. The second term is lagged consumption, whose presence 
can be motivated by habit persistence and/or the presence of rule of thumb consumers (Campbell 
and Mankiw, 1989; Fuhrer, 2000; McCallum and Nelson, 1999). 

“Some degree of arbitrariness is obviously present in the process of recasting the model at a 
different frequency. 

The formulation presented in the text has the advantage of removing most of the seasonality 
in the inflation series, but leads to an inconsistency in the definition of the 1-month real rate 
(which is deflated by inflation expected over the next year, rather than over the next month). We 
used the monthly definition of inflation (i.e. zez In P; — In R-1) and a more theoretically sound 
definition of the real rate (i.e. re Ex [mt+1]) in a previous version of this paper (available online 
from http://www.frbsf.org/economics/conferences/0303/htv.pdf). We ultimately used the former 
specification because the latter has the major disadvantage of leading to a large increase in the 
dimension of the (already quite large) parameter space, since many more lags must be included in 
the system to remove the seasonality of the inflation series. 
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to solve for the rational expectations equilibrium. Since our estimates will include 
also bond prices, we focus on private agents’ perceptions of the monetary policy rule 
followed by the central banks, rather than solving the model under full commitment 
or discretion. Accordingly, the “simple rule” supposedly followed by the central 


bank is to set the nominal short rate according to 


re = (1— p) (8 (Et [mt] — mE) + yee) + pre-a + (3) 


where 7; is the perceived inflation target and 7; is a “monetary policy shock”. 

This is consistent with the formulation in Clarida, Galf and Gertler (1998, 
henceforth CGG), which is a natural benchmark for comparison because the rule 
has been estimated for Germany, the country which we focus on in the empirical 
implementation. The first two terms represent a typical Taylor-type rule (in this 
case forward looking), where the rate responds to deviations of expected inflation 
from the inflation target. The second part of the rule is motivated by interest rate 
smoothing concerns, which seem to be an important empirical feature of the data. 

The main difference with respect to the rule estimated by CGG is that we also 
allow for a time-varying, rather than constant, inflation target 77. We adopt this 
formulation because the Bundesbank modified its “medium term price norm” over 
the sample period used in our analysis and the modifications were public knowledge. 
At the same time, we do not want to impose that the announced price norm was 
“credible,” and reflected in bond prices, by assumption. For this reason, we treat the 
time-varying inflation target 27 as an unobservable variable, which should capture 
markets’ perceptions reflected in equilibrium bond yields. This formulation allows 
us to exploit the full available sample period, without having to assume a break in 
the policy rule at some point in the late seventies, as done by CGG. 

Finally, we need to specify the processes followed by the stochastic variables of 
the model, i.e. the perceived inflation target and the three structural shocks. We 
assume that our 3 macro shocks are serially uncorrelated and normally distributed 


with constant variance. The only factor that we allow to be serially correlated is the 
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unobservable inflation target, which will follow an AR(1) process 


T = bnM-1 + Unt (4) 


where uz is a normal disturbance with constant variance uncorrelated with the 
other structural shocks. 
2.2 A general macroeconomic set-up 


In order to solve the model we write it in the general form 


E 
+ Kr; ep 1,t+1 , (5) 


EX2 t41 X24 0 


H 


Ain = Xit 


where X, is a vector of predetermined variables, Xə includes the variables which 
are not predetermined, re is the policy instrument and E: is a vector of independent, 
normally distributed shocks (see the appendix for the exact definitions of all these 
variables in our example). The short-term rate can be written in the feedback form 
eae | el (6) 

Xat 

This linear structure is nevertheless general enough to accommodate a large 
number of standard macroeconomic models, potentially much more detailed than 
the one we adopt here. The main restriction we impose, for simplicity, is that only 
the short-term interest rate, which is controlled by the central bank, affects the 
macro economy, whereas longer rates do not. 

The solution of the (5)-(6) model can be obtained numerically following stan- 
dard methods. We choose the methodology described in Söderlind (1999), which is 
based on the Schur decomposition. The result are two matrices M and C such that 
Xit Mäi +814 and X94 = CD Consequently, the equilibrium short term 
~ 5The presence of non-predetermined variables in the model implies that there may be multiple 


solutions for some parameter values. We constrain the system to be determinate in the iterative 
process of maximizing the likelihood function. 
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interest rate will be equal tom = A Kit, where A’ = —(Fi+F2C) and F; and 
EF: are partitions of F conformable with Zi: and X2;. Focusing on the short-term 


(policy) interest rate, the solution can be written as 


Te = A'Xi 4 


Xie = Mäi, (7) 


2.3 Adding the term structure to the model 


The system (7) expresses the short term interest rate as a linear function of the 
vector Xı, which in turn follows a first order Gaussian VAR. This structure is 
formally equivalent to that on which affine models are normally built. To derive 
the term structure, we only need to impose the assumption of absence of arbitrage 
opportunities, which guarantees the existence of a risk neutral measure, and to 
specify a process for the stochastic discount factor. 

Behind this formal equivalence, however, our model has the distinguishing fea- 
ture that both the short rate equation and the law of motion of vector Xı have 
been obtained endogenously, as functions of the parameters of the macroeconomic 
model. This contrasts with the standard affine set-up based on unobservable vari- 
ables, where both the short rate equation and the law of motion of the state variables 
are postulated exogenously. 

This feature also differentiates our approach from Ang and Piazzesi’s (2003). 
More specifically, Ang and Piazzesi (2003) still rely on an exogenously postulated 
model of the short-term rate, which they interpret as the monetary policy rule. 
In any macroeconomic model, however, the dynamics of the short term rate will 
be obtained endogenously. We show that this property of macro-models does not 
prevent the specification of a dynamic arbitrage-free term structure model. Provided 
that one’s favorite macroeconomic model can be cast in the linear (5)-(6) form, 
arbitrage-free pricing is possible. 

In fact, rather than building the term structure directly on equations (7), we al- 
low for the possibility to write bond yields as functions of a different vector, Ze, which 


can include any variable in X; or the short term rate. The new vector Z; is defined 
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as Zur = DX;, where D is a selection matrix. Obviously, Z; can also be rewritten as a 
function of the predetermined vector Zu using the result X24 = CX ,. This yields 
Zi = DXi 4, where D is a matrix described in the appendix. Specifically, in the 
empirical application, we choose D so that bond yields are expressed as functions 
of the levels of the macro variables, rather than of their shocks. 

Given the solution equation for the short term interest rate written as a function 
of the Z; vector, re = DZ, we follow the standard dynamic arbitrage-free term 
structure literature and define the (nominal) pricing kernel mz41, which prices all 
nominal bonds in the economy, as mz41 = exp (—rt) e due, where ru is assumed 
to follow the log-normal process un = Y exp (—4A,At — MÉI), 

We then make an assumption on the dynamics of A, the vector of market 
prices of risk associated with the underlying sources of uncertainty in the economy. 
These have commonly been assumed to be constant (in the case of Gaussian models) 
or proportional to the factor volatilities (e.g. Dai and Singleton, 2000), but recent 
research has highlighted the clear benefits in allowing for a more flexible specification 
of the risk prices (e.g. Duffee, 2002; Dai and Singleton, 2002). We therefore assume 


that the market prices of risk are affine in the state vector Z+ 
At An + Au Ze, (8) 


so that the market’s required compensation for bearing risk can vary with the state 
of the economy. 

It should be pointed out here that, in a micro-founded framework, the pricing 
kernel (or stochastic discount factor) would be linked to consumer preferences, rather 
than being postulated exogenously as we do here. The pricing kernel would be 
obtained from the intertemporal consumption Euler equation, essentially consisting 
of the discounted ratios of marginal utility between two consecutive periods, scaled 
by expected inflation in the case of the nominal kernel. In standard consumption- 
based formulations of asset pricing models, the prices of risk would be related to 
the agents’ risk aversion and to the curvature of the indirect utility function with 


respect to the state variables of the problem. We would obtain a micro-founded 
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pricing kernel if we specified a utility function, set A; = 0 and restricted An to be 
consistent with the selected utility function.® 

We prefer our exogenous specification (8) for two main reasons. The first is that 
we want to employ an empirically plausible formulation and the state-dependent 
specification in equation (8) is not straightforward to obtain from first principles.” 
The second reason is that, even if we found a sufficiently flexible formulation of the 
utility function, the yield premia would always be zero in a log-linearized solution 
of the model, such as the one we implicitly adopt here (see also Kim et al., 2003). 
Higher order approximations could obviously be employed to deal with this problem, 
but they would imply leaving the convenient affine world, in which both the bond 
prices and the likelihood can be specified in closed-form. 

In the appendix we show that the reduced form (7) of our macroeconomic 
model, coupled with the aforementioned assumptions on the pricing kernel, implies 
that the continuously compounded yield uff on an n-period zero coupon bond is 
given by 

yy = An + Bi, Zet, (9) 


where the A, and By, matrices can be derived using recursive relations. Stacking 


Consider the example of a standard economy populated by a representative household with 
ry 

utility function over consumption and labour u (Cr, Lt) = EE = L? and production function 

H = A:L?, where A; is a productivity shock such that In A; = pln A-1 + ct. If prices are fully 


flexible, it is easy to show that equilibrium will imply Ct = Y; and, using lower case letters to denote 


the natural log of a variable, yz = ETAC In D + Freya 


The (real) stochastic discount factor Qr t+1 will be given by Qz,441 = Gr = 8 (Yea / Y) 1. 
If we use the definition of the gross interest rate Ry as 1/Ri = Fr [Qr t+1];, take logs, and write 


output in terms of its determinants, we obtain 


1 yo K 2 yo 
REN rt Üe OeVt4+1 
o a Goetz 


where vz41 is a white noise shock with unit variance. 
This formulation is entirely analogous to the one used in the paper, i.e. qtt+1 = —Tt ZALA 
\;£1,t+1. In this simple model, the (constant) prices of risk are given by 


yo 
= o+a(y—-h* 


"Dai (2003) argues that preferences embodying a particular specification of habit formation 
would be consistent with pricing kernel that, to a first order approximation, would be of the form 
(8) with a non-zero An, 
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all yields in a vector Y+, we write the above equations jointly as Y; = A+ B’Z; or, 


equivalently, Y; = A + BX). where Di = BD. 


2.4 Maximum likelihood estimation 


In order to estimate the model, we need to distinguish first between observable and 
unobservable variables in the X; vector. We adopt the approach which is common 
in the finance literature and which involves inverting the relationship between yields 
and unobservable factors (Chen and Scott, 1993). In our case, the method needs 
to be extended to take into account that the observable variables include not just 
the yields, Y+, but also some of the non-predetermined variables. We also use the 
common approach of assuming that some of the yields are imperfectly measured to 
prevent stochastic singularity. 

Using the assumption of orthogonality of measurement error shocks and shocks 
to the unobservable states, we show in the appendix that the log-likelihood function 


to maximize takes the form 


£(0) = -—(T-1) DEE DEIER (27) Su 


r felt 
ier, a)! (BB) Ter 
t=2 i=1 o? 
where A), are the unobservable variables included in the 314 vector, u;” are the 
measurement error shocks, J is a Jacobian matrix defined in the appendix, ©’ is the 
variance-covariance matrix of the four macroeconomic shocks, o: are the standard 
deviations of measurement error shocks, T is the sample size, nm is the number of 
measurement errors and np is the number of variables measured without error. The 
expression for the log-likelihood function above is based on the fact that there is a 
one-to-one mapping from the observable variables (yields, output, and inflation) to 
the unobservables (Xj, and u;”), with J being the Jacobian associated with this 
mapping. 
~ 886 far, we have not imposed any restrictions on the X4¢ vector. In the estimation, however, 


care must be taken to avoid that the unobservable variables included in X1; be linearly dependent. 
If this were the case, the Jacobian matrix would not be invertible. 
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When, as in the model used by Ang and Piazzesi (2003), there is no feedback 
from interest rates to the macro variables, estimation can be performed with a two- 
step procedure. In the more general case analyzed here this is not possible and we 
must estimate the whole system jointly. 

In theory, this is of course preferable. The problem is that the parameter space 
is quite large and therefore the optimization problem of maximizing the likelihood 
function is non-trivial and time consuming. We employ the method of simulated 
annealing, introduced to the econometric literature by Goffe, Ferrier and Rogers 
(1994). The method is developed with an aim towards applications where there 
may be a large number of local optima.? 

One disadvantage of the simulated annealing method is that it does not provide 
us with an estimate of the derivatives, evaluated at the maximum, of the likelihood 
function with respect to the parameter vector, i.e. Oln(£(0)) /00’. These deriv- 
atives are necessary to compute asymptotic estimates of the variance-covariance 
matrix of the parameters. The derivatives could be evaluated numerically, but the 
computation would be based on arbitrarily selected step-lenghts 00, with ensuing 
risks of spurious results because of the highly nonlinear fashion in which the para- 
meters enter the likelihood function. 

To deal with this problem, we rely on analytical results to calculate the Jacobian 
Oln(£ (0)) /00’. The evaluation of the analytical derivatives is quite involved. The 
key steps are described in the appendix. 

-The key parameters of the simulated annealing method were set as follows: Ty = 15; rr = 0.9; 
Nr = 20. The convergence criterion € was set at e = 1.0E — 8. In a preliminary estimation, the 
starting values were taken from CGG’s results (for the policy rule) and from the parameters of an 
unrestricted VAR in output, inflation, and the short term nominal rate. The estimates reported 


in the text correspond to a maximum value of the likelihood function found in a process of 100 
estimations using simulated annealing, starting from randomized initial values. 
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3 An application to German data 


3.1 Data 


Our data set runs from January 1975 to December 1998. The term structure data 
consists of monthly German zero-coupon yields for the maturities 1, 3 and 6 months, 
as well as 1, 3, and 7 years.!9 We assume that the 1-month rate and the 3-year yield 
are perfectly observable, while the other rates are subject to measurement error. 
Yields have been bootstrapped from on an original Bundesbank dataset of end-of- 
month raw prices, coupons and maturities.!! 

Concerning the macro data, we construct the year-on-year inflation series using 
the CPI (all items). For the output gap, we simply follow CGG and detrend the log 
of total industrial production (excluding construction) using a quadratic trend. We 
only deviate from CGG in constructing the series recursively, so that each datapoint 
is obtained by fitting a quadratic trend to the original series up to that point. We 
adopt this approach to ensure that our forecast at time t does not rely on information 
unavailable at that point in time. Both series refer to unified Germany from 1991 
onwards and to West Germany prior to this date. The macroeconomic and term- 


structure series are shown in Figure 1. 


3.2 Estimation results 


To reduce the parameter space in our empirical application, we impose a number of 
restrictions on the coefficients of the market prices of risk. In the general set-up, we 
showed that the risk prices can be specified as Au = An + Au Z. In our application, 
Z+ includes the perceived inflation target and contemporaneous and lagged values of 
inflation, output and the short term rate. Given Zu, Au can obviously have nonzero 
elements only corresponding to time t variables, as lagged variables are no longer 
subject to surprise changes. This leaves only four potentially non-zero rows in the An 
and A, matrices, corresponding to the perceived inflation target, the policy interest 
~ 10We do not use 10-year bonds because these are only available without breaks as of April 1986. 

Che methodology is equivalent to that employed by Fama and Bliss (1987). We wish to thank 


Thomas Werner for providing us with the raw data and Vincent Brousseau for bootstrapping the 
term structures of zero-coupon rates. 
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rate, inflation and the output gap. Next, we restrict A, in the sense of allowing 
interactions only between prices of risk of contemporaneous variables, which leaves 
us with a 4 x 4 non-zero submatrix in Au. Finally, we follow Duffee (2002) and 
set to zero all entries whose elements have a t-statistic lower than 1 in preliminary 
estimations.'? As a result, we are left with the following non-zero elements in the 


matrices of prices of risk 


Aol 0 0 A13 À14 Ti 

Ao2 A21 A22 A23 U rt 
Se $ 

A Aan A32 Aen 0 Tt 

A 0 Ag 0 Au Lt 


3.2.1 Parameter estimates 


Table 1 presents the parameter estimates with associated asymptotic standard errors 
(based on the analytical outer-product estimate of the information matrix). 

The results are broadly consistent with the evidence of Clarida, Galf and 
Gertler (1998) regarding the Taylor rule in Germany and, as far as the other macro- 
parameters are concerned, with existing evidence based on structural models or 
identified VARs. 

For example, our point estimate of the degree of forward-lookingness of infla- 
tion (4) is within the range of values found by Jondeau and Le Bihan (2001), who 
estimate on German data a Phillips curve based on quarterly data using a variety of 
specifications and two different estimation methods. Kremer, Lombardo and Werner 
(2003), who estimate a structural macroeconomic model with explicit microfounda- 
tions, estimate a much higher value of u„. Their estimate, however, is not directly 
comparable to ours due to the fact that they capture the persistence of inflation 
through highly persistent exogenous shocks (whereas our shocks are white noise). 

A result which casts doubts on the ability of our macro-model to provide an 
~ 12These preliminary estimations involved first estimating the model with full 4 x 1 Ae and 4 x 4 
An matrices. Then, the A parameter with lowest t-statistic (below 1) was fixed at zero. The model 


was reestimated, and again the A parameter with lowest t-statistic was constrained to be zero. This 
process was repeated until no A parameters with t-statistics below 1 remained. 
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accurate characterization of the dynamics of output and inflation in Germany is 
that the elasticity of inflation to the output gap is very small (ôs = 0.0004 and 
insignificantly different from zero). This is not entirely surprising. Jondeau and Le 
Bihan (2001) also find values of 6, close to zero for some specification/estimation 
method (Kremer, Lombardo and Werner, 2003, calibrate, rather than estimate, this 
parameter). Identified VARs estimated at the monthly frequency (e.g. Sims, 1992) 
also tend to find a very small and insignificant responses of inflation to, e.g., mone- 
tary policy shocks, which is consistent with our results of a very small 6, and also 
a small ¢,.. 

To assess whether our macro-parameter estimates are affected by our inclusion 
of term structure information in the model, we re-estimated the macroeconomic 
model separately. In order to work with a more conventional set-up, we also elim- 
inated the stochastic inflation target from the policy rule and replaced it with the 
Bundesbank’s announced price norm. Apart from a small increase in ¢, from 0.03 
to 0.06, the other parameter estimates (including ôy) were virtually unchanged. 

The macro-model performance may be affected by the fact that volatile, monthly 
data are noisy and make it harder to identify the link between inflation, output and 
interest rates. Another possibility is that our output gap definition, which plays a 
crucial role in the analysis, is an imperfect proxy for the theoretical notion of real 
marginal costs. Or else, as already emphasized, our 2-variable macro-model may be 
too parsimonious to describe German macroeconomic dynamics, which are possibly 
affected also by variables such as the exchange rate or, as in Kremer, Lombardo 
and Werner (2003), a monetary aggregate. Since our main interest is not that of 
finding the macroeconomic model most suited for German policy analysis, we do not 
perform further specification search. We only test for a potential missing variable 
bias by examining the residuals’ autocorrelation. We find little evidence of serial 
correlation in our preferred specification. 1? 

As to the other parameters, the autocorrelation coefficient of the inflation target 


'3More precisely, looking at the correlograms of the estimated residuals we find no evidence of 
statistically significant first or higher order correlation in the output and inflation equations. 
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process is very close to 1.14 Concerning the term structure, our estimates of the 
standard deviations of the measurement errors are between 23 basis points for the 
3-month rate and 28 basis points for the other yields. These values are broadly in 
line with the results of models based solely on unobservable factors and also those of 
an unrestricted VAR including inflation, the output gap and the bond yields.! The 
standard errors of the 1-month and 3-month rate equations of the VAR are equal 
to 43 and 32 basis points, respectively, compared to 48 and 23 in our model; for 1- 
year and 7-year yields, the VAR equations have a standard error of 29 and 24 basis 
points, respectively, compared to 28 and 28 in our model. Obviously, our model 
has the advantage of describing, at the same time, the yields for all other possible 
maturities (and it also does better than the VAR at fitting output and inflation). 
Finally, one of the benefits of our model is that of providing us with a measure 

of the central bank’s inflation target as reflected in the prices of long term bonds. 
One of the tests of the model is therefore to check whether the filtered series “looks” 
reasonable. For this purpose, Figure 2 compares it to the Bundesbank medium term 
price norm.!® The two series are quite close to each other in the volatile seventies 
and in the sharp decline of the beginning of the eighties. A large discrepancy can be 
observed mostly at the beginning of the nineties, when the estimated target increases 
sharply while the price norm remains unchanged. The increase in the estimated tar- 
get is, however, not unreasonable, as it coincides with an increase in actual inflation 
following the expansionary policies that accompanied German unification.!’ At the 
same time, we cannot exclude the possibility that the variability of the target may be 
overestimated in order to induce inflation persistence present in the data, which our 
stylized model may be unable to capture endogenously. Nevertheless, the perceived 
inflation target is less variable than actual inflation, both in terms of its sample 
standard deviation and of its minimum and maximum sample values. 

14This parameter is constrained to be strictly smaller than 1 in the estimation. 

15 The VAR is estimated over the same sample period and includes 3 lags of the variables. 
Until 1981 and from 1997 to 1998, the Bundesbank actually announced a range, rather a point 
value, for the price norm. In these years, the mid-point of the range is displayed in Figure 2. No 
values were announced pre-1976 and in 1979. 

'7Tn spite of the unchanged price norm, this may have sparked fears of a waning in the Bundesbank 
anti-inflationary determination because of domestic — due to unification — and European-wide — due 


to the impact of any monetary policy tightening on ERM partner countries — political pressures (see 
Issing, 2003, for a concise account of the Bundesbank’s policy at the time of German unification). 
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3.2.2 Impulse response functions 


Our structural model allows us to compute impulse response functions of macro 
variables and yields to the underlying macro shocks. 

Figures 3 to 6 show the impulse responses of selected variables to the structural 
shocks. The responses of the macroeconomic variables and of the short term interest 
rate are broadly in line with existing VAR evidence based on German (monthly) 
dataand we will not delve on them here. We concentrate instead on the responses 
of yields. 

We start from Figure 3, which displays the impulse responses to a shock to the 
perceived inflation target, which increases on impact by approximately 0.2 percent- 
age points. The shock is obviously expansionary and very persistent, due to the high 
serial correlation of the inflation target process. The response of the yield curve is 
an almost parallel and very persistent upward shift at all maturities, except the very 
short ones (which move slowly because of the high interest rate smoothing coefficient 
in the policy rule). The size of the shift corresponds roughly to that of the initial 
inflation target shock and it is significantly different from zero for maturities around 
1-year. 

Figure 4 shows the effect of a 45 basis points increase in the 1-month interest 
rate because of a monetary policy shock (the disturbance 7,). The response of the 
yield curve is decreasing in the maturity of yields, which factor in the slow return 
to baseline of the policy rate. Hence, a monetary policy shock tends to cause a 
statistically significant change in the slope of the yield curve. The shape of this 
response is qualitatively similar to that obtained by Evans and Marshall (1996) for 
the US. 

An inflation shock, shown in Figures 5, tends to increase the curvature of the 
yield curve. Yields move little and slowly at the short end, more around the 1- 
year maturity, then little again at the long end. While statistically significant for 
maturities below 7 years, the responses appear to be very small from a quantitative 


viewpoint. 
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Finally, Figures 6 shows the impulse responses to an output shock. Due to 
the small policy response, the yield curve increases little, but significantly, over 
maturities up to 1 year. Yields on 3 and 7-year bonds, however, fall as a result of 
the shock and in spite of the fact that the response of the short-term rate always 
remains above the baseline. This surprising pattern is to a large extent shaped by 


the dynamics of risk premia. 


3.3 Macro shocks and risk premia 


Another advantage of our joint treatment of macroeconomics and term-structure 
dynamics is that we are able to derive the impulse response of theoretical risk premia 
to macro shocks, including the monetary policy shock. These are shown in Figure 
T. 

The inflation target shock is immediately followed by an increase of the yield 
premium for maturities up to 4 years, with a peak effect of 10 basis points at the 
l-year maturity. The premium then turns negative for longer maturities. Such 
increase in the yield premium is highly significant from an economic viewpoint, as it 
plays a large quantitative role in shaping the total yield response displayed in Figure 
3. 

The monetary policy shock gives rise to a large fall, on impact, at the short 
end of the term structure of yield premia, thus reducing significantly the size of 
the impact response of the yields. The impact response of the 1-year yield to the 
monetary policy shock, for example, would increase by one half if yield premia were 
set equal to a constant. 

Similar considerations hold for the impact response of yield premia to inflation 
and output shocks. The latter is notable, since the premia embody most of the 
action in the response. The impact response of the 7-year rate, for example, would 
change sign and essentially maintain the same in absolute value, if risk premia were 
constant. 

We conclude that, in general, the dynamics of yield premia have a nonnegligible 
effect on the impulse responses of yields to all macroeconomic shocks. An interpre- 
tation of the yield responses based on the expectations hypothesis may therefore be 


significantly biased. 
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The general features of the yield premia are that their level and volatility are 
increasing in maturity. The premia also tend to be decreasing over the sample in 
parallel to the fall in inflation, but then shoot up again, temporarily, in 1992-93. To 
investigate their determinants more closely (using equation (15) in the appendix), we 
can decompose the premia in the components due to risk of changes in the inflation 
target, in the short-term rate, in inflation and in the output gap.!8 We focus here 
on the most important components for 1- and 7-year maturities. 

The most striking outcome of this decomposition is that premia linked to in- 
flation risk are almost perfectly constant over time and negligible in size across 
maturities. Even at their peaks, they never reach the level of 10 basis points. This 
number should be compared, for example, to the maximum level of 1 percentage 
point reached by the premium due to output gap uncertainty for 7-year bonds. 

Variations in yield premia arise by and large from fluctuations in the other three 
variables, with an importance that changes across maturities. Figure 8 shows that 
at the 1-year horizon, the largest fraction of the time-varying yield premium is due 
to interest rate risk, i.e. the possibility of monetary policy surprises. The price of 
interest rate risk, in turn, is decreasing in the level of the interest rate: when the 
latter is very high, yield premia are lower than average and 1-year bonds appear to be 
a very appealing form of investment; when interest rates are low, on the contrary, the 
risk of unexpected changes in the short-term rate has a higher price and 1-year bond 
command a higher than average premium. The second most important component 
of the time varying yield premium at 1-year maturities is inflation target risk. The 
target premium is increasing in the level of the inflation target. A high inflation 
target makes 1-year bonds riskier and increases the premium investors require to 
hold them. 

At the long 7-year horizon, Figure 8 shows that the time varying component 
of the yield premium is almost entirely due to inflation target risk until the end of 
~ 18Thig decomposition is not exact, because the term premium is also affected by the lags of 
inflation, output and the interest rate. We disregard these additional effects for two reasons. First, 


given our assumption on the prices of risk A, they are due to convexity effects, rather than a pure 
risk premium. Second, they are quantitatively minor. 
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1988. At this maturity, the inflation target premium is negatively correlated with 
the level of the inflation target. When the target is high, the yield premium is lower 
than average and investors are relatively more willing to hold 7-year bonds. This 
may be taken as a signal of investors’ confidence in the ultimate return to a low 
inflation target environment. As of 1989, with inflation and the policy interest rate 
increasing after the German unification and the recession of 1992-93 ensuing, the 
variable yield premium becomes significantly affected also by output gap risk. In 
other words, booms tend to make investors more willing to hold long term bonds, 


while investors require a larger bond premium during recessions. 


4 Forecasting 


The forecasting performance is a particularly interesting test of our macroeconomic- 
based term-structure model. Due to the relatively large number of parameters that 
needs to be estimated, the model could be expected to perform poorly with respect to 
more parsimonious representations of the data. In fact, the random walk model has 
been shown to provide yield forecasts that are particularly difficult to beat (Duffee, 
2002). An important test of our model is therefore to check whether the information 
contained in macro variables can improve the performance of a standard essentially- 
affine model including only term-structure information. For completeness, we also 
check whether the inclusion of yields in the information set can improve the perfor- 
mance of the macro-only model in terms of forecasting the macro variables. 

The forecasting tests for macroeconomic variables and yields are presented in 
turn in the next two sections. Our results suggest that term structure information 
helps little in forecasting macroeconomic variables. Our structural framework in- 
cluding macroeconomic variables does, however, help to forecast yields. The out-of- 
sample forecasting performance of our model up to 12-month ahead is almost always 
superior to all the alternatives we consider, and the difference is often statistically 


significant. 
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4.1 Do yields help to forecast macroeconomic variables? 


Given the imperfect ability of our macroeconomic model to describe the joint dynam- 
ics of German macroeconomic variables, we do not expect it to be very successful in 
forecasting inflation and the output gap. This is consistent with existing evidence. 
In a thorough study of inflation forecasting in the G7 countries, for example, Canova 
(2002) concludes that theory-based models are not always better than atheoretical 
univariate models. 

Our test on macroeconomic variables is therefore very focused on assessing 
whether including yields in the analysis can help in forecasting. We evaluate the 
out-of-sample forecasting performance of our macro-term structure model (HTV) 
relative to the macro-only model (HTV-M) used in Section 3.2.1. Both models are 
estimated over the period February 1975 - December 1994, and output gap and 
inflation forecasts up to 12 months ahead are produced for each month over the 
period January 1995 to December 1998. We update the information set each month, 
but we do not reestimate the model. Instead, we rely on the estimates up until end- 
1994. We choose this approach to limit the computational burden of the exercise 
(this applies in particular for the yield forecasts in the next section). All results are 
therefore based on the same estimated parameters. 

The root mean squared errors (RMSEs) are presented in Table 2. Lower val- 
ues of the RMSE denote better forecasts, and the best forecast at each horizon 
is highlighted in bold. The table shows that the full model including term struc- 
ture information and the stochastic inflation target does marginally better than the 
macro-only model at forecasting inflation. The latter model, however, prevails as 
far as output forecasts are concerned. 

We conclude that yields are unlikely to provide useful information for macroeco- 
nomic forecasting within our framework. This result may be due to our assumption 


that long term yields do not affect the dynamics of inflation and the output gap. 
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4.2 Do macroeconomic variables help to forecast yields? 


To assess the yields forecasting performance of our model, we compare it to a number 
of benchmarks. 

The first is the random walk. In addition, we also consider forecasts based 
on three other models. One is a canonical An (3) essentially affine model based on 
unobservable factors.! Provided that risk premia are specified to be linear functions 
of the states, Duffee (2002) finds this model most successful in the class of admissible 
affine three factor models in terms of forecasting US yields. Apart from providing 
a benchmark for comparison, our results on the Ao (3) model are of independent 
interest, since they highlight the performance of this model on a different data-set. 
The second model we take into account is the Ang and Piazzesi (2003) model, which 
we reestimate on our data-set. Based on Ang and Piazzesi’s results, we use their 
favorite “Macro model” in this exercise, i.e. a model in which the interest rate 
responds to current inflation and output gap, as well as to 3 unobservable factors. A 
potentially important difference in our application of their model, however, is that we 
use inflation and the output gap directly in the estimation, rather than the principal 
components of real and nominal variables employed by Ang and Piazzesi (2003), 
thereby facilitating comparison to our results. Finally, we use an unrestricted VAR 
including all the variables in our structural model, in order to gauge the importance 
of structural and no-arbitrage restrictions to improve the performance of our model. 

As in the previous section, the out-of-sample forecasting performance of the 
models is based on estimates up to end-1994, and 1- to Listen ahead forecasts 
are produced for all yields used in the estimation over the period January 1995 to 
December 1998. The RMSEs of this forecast evaluation exercise are summarized 
in Tables 3, with the best forecast at each maturity/horizon highlighted in bold. 
The exercise shows that our model performs better than the alternatives for all 
maturities, at least beyond the very shortest forecast horizon. In particular, our 
model beats the predictions of the random walk benchmark in almost all cases. Table 


4, which displays the trace MSE statistic — a multivariate summary measure of the 


'°For a definition of the Ao (3) class of affine models, see Dai and Singleton (2000). 
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forecasting performance across yields for each horizon — confirms this picture 2 


To understand the reasons for this success, compare first the performance of the 
Ao (3) model in Table 3 to that of the VAR. The former model includes no-arbitrage 
restrictions and, as a result, it appears to be more efficient at forecasting long yields, 
especially at longer forecasting horizon. The An (3) model, however, is not always 
superior to the VAR, which is a first suggestion that macroeconomic information 
could be important in forecasting yields. The AP may be expected to improve the 
performance of the Ao (3) model, because it includes macroeconomic information on 
top of the no-arbitrage restrictions. The AP model includes, however, a very large 
number of parameters to estimate, since it is based on a reduced-form representation 
of the macroeconomic variables. This may be the cause for its less satisfactory 
performance over forecasting horizons beyond 1 month. Its good performance in 
1-step ahead forecasts is, incidentally, consistent with the results reported by Ang 
and Piazzesi (2003). Our model appears to strike a good balance in incorporating 
macroeconomic information without becoming overparameterized. 

Concerning, more specifically, the market prices of risk, a crucial role in affecting 
the forecasting performance of our model is played by risk premia associated to 
inflation target risk. If we re-estimate our model restricting to zero the elements 
in equation (8) associated to the inflation target, i.e. Agi and Aan, the forecasting 
performance of the model worsens dramatically, especially for long maturities. This 
appears to be consistent with the evidence on the main components of the risk 
premia presented in section 3.3. 

In order to formally test the out-of-sample yield forecasting performance of our 
model, we apply White’s (2000) “reality check” test. This test, which builds on the 
work of Diebold and Mariano (1995) and West (1996), involves examining whether 
the expected value of the forecast loss (e.g. the squared forecast error) of one or 
several models is significantly greater than the forecast loss of a benchmark model. 
We choose this test mainly because, in contrast to many other forecast performance 
~ 20-The trace MSE statistic has been used, for instance, by Christoffersen and Diebold (1998). For 
each forecast horizon, it is simply computed as the trace of the covariance matrix of the forecast 


errors of all yields considered. Hence, a lower trace MSE statistic signals more accurate forecasts 
across yields. 
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tests, it tests for superior predictive ability rather than equal predictive ability. 
Moreover, White’s test allows us to examine whether a specific model is significantly 
outperformed by any model among a number of alternatives, whereas other tests 
typically do not permit this. 

We implement the test in two ways. First, we take our model as the benchmark 
and ask whether any of the four alternative models is able to produce forecasts that 
are significantly superior to our model across each of the five maturities and 12 
forecast horizons considered. We found that in only 3 out of the 60 cases could we 
reject the null hypothesis that none of the four models is better than our model.?! 
While encouraging, this result does not necessarily imply that our model is superior 
to the alternatives. To test this, we instead turn around the null hypothesis and 
proceed to test for superior predictive ability of our model vis-a-vis each of the 
four alternative models separately. In other words, we implement four pairwise 
comparisons, where each of the four alternative models are used as the benchmark 
model against the HTV model. The results are displayed in Table 5, where bold 
figures indicate rejection of the null that our model does not have superior predictive 
ability compared to the benchmark used, at the 5% level. In over 60% of the cases we 
reject the null, meaning that for most of the combinations of maturities and forecast 
horizons considered here, the forecasting performance of our model is significantly 
better than the alternatives. Looking at the results in more detail, we see that, 
somewhat surprisingly, the VAR model seems to be harder to beat than the other 
alternatives, although for longer horizons and maturities the HTV model consistently 
outperforms the VAR. With respect to the performance of our model at different 
forecast horizons, we seem to do roughly equally well across all horizons, except for 
the one-month ahead case, where the null is rejected less often. 

We therefore conclude that the inclusion of macroeconomic variables within 
a structural framework contributes to sharpening our ability of forecasting yields 
accurately out of sample. The improvement is due both to the inclusion of addi- 


tional information in the model, and to the structural restrictions imposed on its 


macroeconomic and term structure sections. 


21 These results are not reported, but are available on request. 
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5 Expectations hypothesis tests 


According to the expectations hypothesis, the yield on an n-period zero-coupon 
bond should increase when the spread between the same yield and the short term 
rate (the “slope of the yield curve”) widens. In fact, the projection of the yield 
change aE — y? on the yield spread (y? — r+) /(n — 1) should yield a coefficient 
of 1. A number of empirical tests of this implication of the theory have, however, 
found a negative relationship. This pattern represents a puzzle for the expectations 
hypothesis, and it appears to be particularly clear for United States data. The 
relevant regression coefficient can be as big as —5 for 10-year bonds, according to 
e.g. Campbell and Shiller (1991), while the expectations hypothesis predicts a value 
of one for all maturities. 

One interpretation of these results is that the large deviation from 1 in the 
estimated coefficient on the yield spread is due to large and time varying risk pre- 
mia (not permitted by the expectations hypothesis). Using a highly stylized model, 
McCallum (1994) conjectures that an exogenous, stochastic term premium is, in 
principle, capable of causing deviations from 1 in the slope coefficient of the afore- 
mentioned regression. The actual size of the deviation will depend on both the 
stochastic properties of the term premium (see also Roberds and Whiteman, 1999) 
and the monetary policy rule followed by the central bank. These papers, however, 
work by “reverse engineering.” Given the results of projections of the yield change on 
the yield spread, they derive the properties that risk premia should have to explain 
those results. This is different from deriving the risk premia from a certain model 
and checking ex-post if they are capable of solving the expectations puzzle. 

In this section, we follow the latter strategy. We do not test if the yield premia 
consistent with our model are capable of solving the expectations puzzle for some 
parameter values. This is likely to be the case, given that our model includes a 
relatively flexible specification of the market prices of risk. We ask instead a more 


stringent question, namely whether the premia generated by our model can solve 
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the expectations puzzle for the specific set of parameter values which maximizes the 
likelihood. 

In so doing, we follow closely Dai and Singleton (2002) who ask the same ques- 
tion within a number of dynamic affine term structure models based on unobservable 
factors. More specifically, we ask whether the model-implied, population coefficients 


o, in the regression 
thor — y = const. + o, (yz — rt) / (n — 1) + residual (10) 


match the values obtained from an OLS regression on actual yield data. The pop- 
ulation coefficient are obtained assuming that the model parameters are true and 
then deriving the o. coefficients analytically based on the stochastic properties of 
the model.?? Following Dai and Singleton (2002), we also examine the small-sample 
counterparts of these coefficients. Some correction for small sample bias is desirable 
because of the persistent nature of yields. For this purpose, we generate 1000 sam- 
ples of the same length of our data (287) and calculate the mean estimate of the d, 
coefficients. 

Dai and Singleton (2002) denote the above test as LPY (i). LPY (i) is a test of 
the capacity of the model to replicate the historical dynamics of yields as generated 
by a combination of the dynamics of risk premia and expectations of future short 
rates. As already emphasized, a successful model should be capable of generating 
the negative intercept coefficient of Campbell and Shiller-type regressions. 

In addition to LPY (i), Dai and Singleton (2002) also suggest running a sec- 
ond sort of test, defined as LPY (ii), which focuses on the realism of the dynamic 
properties of risk premia. If the model captures these dynamics well, a Campbell 
and Shiller-type regression based on risk-premium-adjusted yield changes should re- 


cover the coefficient of unity consistent with the expectations hypothesis. LPY (ii) 


cov(y a", (vt —re)/(m—1)) 


*? This amounts to evaluating analytically d. = een (CN eT ; 
Le = 
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therefore tests that the sample coefficient ø% in the regression 
ea — yf +enst/ (n — 1) = const. + oF, (yg — re) / (n — 1) + residual (11) 


is equal to its population value of 1 (in the above regression, en,t is the excess holding 
period return én, = Er [ln (ni. (nl — rs] — see appendix). 

Dai and Singleton (2002) show that an affine 3-factor model with Gaussian in- 
novations and including a risk-premium specification of the type suggested by Duffee 
(2002) scores extremely well in terms of both LPY (i) and LPY (ii). Our model also 
includes a flexible specification of the risk-premium as in Duffee (2002). Unlike in 
pure finance models, however, our risk-premia are partly functions of observable vari- 
ables, namely output and inflation. This feature represents an additional constraint, 


which makes the LPY tests more stringent than in the pure finance literature. 


5.1 LPY(i) 


Since the evidence on Campbell and Shiller-type regressions based on European data 
is less compelling than for the US (e.g. Hardouvelis, 1994, Gerlach and Smets, 1997, 
Bekaert and Hodrick, 2001), we start by replicating Campbell and Shiller’s analysis 
on our data. The results of the sample estimates of equation (10) are shown in 
Figures 9 and 10 as dots under the label “Sample”. Consistently with the puzzle, 
the estimated intercept coefficient is always negative and decreasing in maturity. We 
confirm, however, that the puzzle appears less severe for German yield data: the 
estimated coefficient hovers around —0.7 for 7-year yields, compared to a value of 
less than —3 reported by Dai and Singleton (2002) for US 7-year yields. 

In Figure 9 we show the results of the LPY (i) test. The population coefficients 
follow quite closely the pattern of the sample coefficients, although less so for short 
maturities. The population coefficients also have the downward sloping feature em- 
phasized also by Dai and Singleton. The small-sample values of the d, coefficients 
(labelled “Model-implied MC” in Figure 9 and drawn together with 95% confidence 
bands of their small-sample distribution) confirm and strengthen this result. Our 
model appears to match strikingly well the pattern of the sample coefficients for 


essentially all maturities included in the regression. 
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The success of the model in matching LPY (i) depends crucially on our assump- 
tions related to the market prices of risk. Our parameterization of the An matrix 
permits variations of the prices attached to the various sources of risk depending 
on the level of the state variables of the model. For example, the risk premium 
required for the possible occurrence of inflation target shocks varies with the levels 
of inflation and the output gap (see the first row of the Au matrix). In fact, it turns 
out that the statistically significant dimension of the inflation target premium is not 
related to the occurrence of “own shocks” (the first element in the matrix is zero). 
What matters is whether inflation and the output gap are high when the target is 
also high because of past inflation target shocks. 

In specifications not allowing for such interactions — for example if the A; matrix 
were diagonal — we experienced a worsening of the the performance of our model 
in terms of the LPY tests. The importance of the interactions generated by the 
off-diagonal terms in the Au matrix is related to the fact that these increase the 
persistence of the yield premia. Once the premium related to inflation target shocks 
has gone up, it will possibly remain high not only because of the persistence of the 
inflation target, but also because of increases in the output gap or inflation driven 
by any other shock in the system. The persistence of yield premia, in turn, is crucial 
to generate significant deviations in the yields levels from the values consistent with 


the expectations hypothesis. 


5.2 LPY (ii) 


Figure 10 shows the results of the LPY (ii) tests. Once again, the model does 
remarkably well in fitting the data. The risk-premium correction always goes in 
the right direction and the model can generate a coefficient very close to unity for 
maturities of 4 years or longer. 

For shorter maturities the model does less well, but we still recover coefficients 
that are positive and larger than 0.5, which is a dramatic improvement with respect 
to the implications of the expectations hypothesis. The reduced degree of success of 


the model at the short end of the yield curve is also consistent with standard results 
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that 3-factor models are unable to capture short-lived money-market dynamics and 
that a fourth factor may be necessary for this purpose. Alternatively, such dynamics 
may be captured allowing for jumps in the short-term rate, as in Piazzesi (2003). 
To summarize, our model appears to do as well as the essentially affine Ao (3) 
class in tests of the expectations hypothesis, in spite of the further constraints im- 
posed by the dependence of risk premia on observed variables. The results of LPY (i) 
are very positive, in that the model can replicate the estimated coefficient of Camp- 
bell and Shiller-type regressions at all maturities. The test of LPY (i) are also 


positive, and especially so for long maturities. 


6 Conclusions 


This paper presents a general set-up allowing to jointly model and estimate a 
macroeconomic-plus-term-structure model. The model extends the term-structure 
literature, since it shows how to derive bond prices using no-arbitrage conditions 
based on an explicit structural macroeconomic model, including both forward-looking 
and backward-looking elements. At the same time, we extend the macroeconomic 
literature by studying the term structure implications of a standard macro model 
within a dynamic no-arbitrage framework. 

In an empirical application, we show that there are synergies to be exploited 
from current advances in macroeconomic and term-structure modelling. The two 
approaches can be seen as complementary and, when used jointly, give rise to sensible 
results. Notably, we show that our estimates of macroeconomic parameters, that are 
partly determined by the term-structure data, are consistent with those that would 
be estimated using only macroeconomic information. At the same time, our model’s 
explanatory power for the term-structure is comparable to that of term-structure 
models based only on unobservable variables. 

We assess the performance of our model mainly along two dimensions: fore- 
casting and ability to solve the expectations hypothesis puzzle. 

While yields do not seem to provide useful additional information in forecasting 


macroeconomic variables, our model performs very well in forecasting yields. We 
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argue that this is due to both the inclusion of macroeconomic variables in the infor- 
mation set and to the imposition of a large number of no-arbitrage and structural 
restrictions. 

Our macro-based term-structure model can also match features of yield curve 
data which represent a puzzle for the expectations hypothesis. These results confirm 
that the dynamics of stochastic risk premia are important determinants of yield 
dynamics, and that all such dynamics can be ultimately reconducted to underlying 


macroeconomic dynamics within a consistent framework. 
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A Appendix 


A.1 State-space form 


We write the model (1)-(3) in the state-space form (5)-(6), we define the vectors 
Xit, Kar and F and the matrix H as follows: 


* gr 20 1 
Xit Lo, 1, Tt—2, ETA AE Tt-2, Tt ae Ms Er Et Tt—1] j 


/ 
Kor = Lëkëeun, — e141, Lt, Baton: Rote, Te] , 


/ 
X, = E si, 


F 


TEEN 


1x11 
Hy = 


0 I 0 
2x3 2x2 2x6 


de Ga ia 


0 
4x11 


O 1 0 
1x11 1x12 
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= 0 
My T 1x6 Hx 
(0) 
Hə, = 11x11 
21 12(1— u) ep 12 g 
0 6, 0 -= 
1x3 Hr 1x2 Ur 1x2 
0 
11x11 
1 12 12¢, 0 
1x11 He He 1x11 
oa jie 
x x 
Hə = D 
0 126% 12 
1x11 Hr 1x11 #r 
0 I 0 
11x12 11x11 11x1 


where Ce = [nena and ôr = äi, Au, O73] ; i 
g 12¢,. — 1 
Finally, we define K = | e i a | and Eu = Kä , 
/ 
where €) 441 = D Un*t41) Matz), MÉ ns Uet41,0| . We can therefore write the 
j X b i 


system as 


Xi = QX; + éni (12) 


where Q = H — KF. 


A.2 Bond prices 


For the pricing of bonds, we work with the transformed vector Z+ defined as Z; = 


* / ` : 
Lo, 1, Tt—2, Lt—3, TWt-1, Mt—-2, Tt 3 TE, Ty EN . Using the solution Xə = CX, 


Z can be written as Zu = DXi, where D is 


(tea) 
TXT 7x4 


-FG 


/ 
G= Du C | and Or: denotes row j of the matrix C. 
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Given the definition of r; and sp1, the pricing kernel me) = exp (—r;) oa 


can be written as 


1 
JAA Nena) 5 (13) 


Mt+1 = exp ( A3 


where we used r; = A Z; with A’ = | 0,1, 0 | 
1x7 1x3 
We know that this set-up will deliver bond prices that are exponential affine 
functions of X14. Since Ze is an affine transformation of X1, we can write the bond 
prices as 


pt = exp (An + Bi, Ze) (14) 


where the coefficients A, and Bn have to be determined. 

Note first that the price of a one-period bond at time t is p} = Ey [m1] = 
exp (A) Zt), so that A; = 0 and Bı = —A. We can now use the pricing kernel 
(13) and the postulated form of bond prices (14) to rewrite the equation for the 


price of an (n + 1)-period bond p?! = E [mpk] as 


7 T EP ae ee eee eee eee 
po Ses (4, - B DEM + 5B, DI'N Bn + (B, ÔMÔ "sA e BDS) d 
where we also used the properties of a lognormal variable zue such that E [exp (a + be1141)] = 
exp (a+ 3b? var [61t41]). The bond-pricing coefficients for any maturity n can there- 


fore be found using the recursion 


_ ners T ow oe 
Au = An- B DEM + 5B,DID'D'B,, 


1 =z! 


Al A. 


D, = BD (Mb 
initialized at Au = 0 and Bı = —A. 


A.3 Likelihood function 


To implement ML estimation of the model, we first partition the state vector X1, 
into a vector X{, that includes only unobservable variables and a vector X9, of 


observable variables. Similarly, we define a vector X5, of observables from Xə. 
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Moreover, to prevent stochastic singularity, we assume that some of the yields are 
subject to measurement errors, that are assumed to be serially uncorrelated and 
mean-zero. We denote these by Y7”, while Y? will denote the (perfectly observed) 
remaining yields. If we denote by np the number of unobservable variables (i.e. the 
dimension of vector X{,), then the sum of the dimensions of X$ , and Y? must also 
equal np. Correspondingly, we will denote by nm the number of variables subject to 
measurement error (i.e. the dimension of NI"), 


In our application, we have 


lt = |Tt—1, Tt—2, Lt—3, Mt-1, Tt—2, START? il’, 
Lit = T$, Np ER EFI, 

2t = Li Tl, 

v= ma] 

Y ea [on a | 


where y;’ denotes the yield on a zero-coupon bond with n-month maturity. 

Next, we follow Chen and Scott (1993), Duffee (2002), and Ang and Piazzesi 
(2003), among others, and use the perfectly observed yields and macro variables 
to back out the vector of unobservable state variables, Xï. To do this, we use 
the fact that Y? can be expressed as Y? = A? + B?X?, + BPX}, where the 
superscript p denotes the selection of factor loadings corresponding to Y”. Similarly, 

/ 
C 


given the relationship X5, = CX where C = C (2 EIE we can write 


2% = C°x?, ar "XY. 


/ 
{12,.} 


Given the vector of parameters 0, these equations can be inverted to form an 
implied vector Xt. More specifically, let W+, denote the vector stacking Y? and 
/ 
Si ie. We=[-y? xy, |- We obtain 


i l AP Be 
tt =J W= a = lt 
0 C° 
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where the Jacobian term is given by 


Finally, given the vector x. implied yields for the remaining nm bonds can be 
computed using Wi = A™ + Bex, + Burx Te Where the superscript m denotes 
the selection of factor loadings corresponding to Y™. In general, these implied yields 
will not exactly correspond to the observed yields. The difference produces the 
vector of measurement errors, u?” = NI — WI, which is assumed to have a constant 
diagonal variance covariance matrix with element 7 given by Ge 

To compute the log-likelihood value, we start from the knowledge that the 
one-period ahead conditional distribution of the unobservable state variables has 
a multivariate normal distribution fx» (Xt, | Xý) This distribution is known, 
since the conditional mean of A7 is given by the theoretical model and its variance- 
covariance matrix XY’ is assumed to be constant and diagonal. The distribution of 


W, conditional on lagged values of the observable states, is then?’ 


o 1 çu u 
fw (W: | Kg = ext (Xi, | Risa) : 


Assuming that the yield measurement errors are jointly normal with distribu- 
tion fum (u7) , the log-likelihood of observation t will be given by the sum In fw (W; | X?,) + 
In fum (u), which can be written as 


T 1)», 


7 ( Pa , 
> ëmgi = -(F-1)njJ|- (27) - —— n |E=" 


2 
n T n m 
(T — 1)nm T-1 3 a 1 D m (ur) 


23 Note that EI, contains only lagged values of the output gap, inflation, and the short-term 
interest rate. 
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Our maximum likelihood estimate is the vector 0* which maximizes the above 


expression. 


A.4 Analytical derivatives 


The calculation of the analytical derivatives of the log-likelihood function with re- 
spect to the parameter vector involves two key steps. First, the derivatives of the 
A and B matrices with respect to the M and C matrices; second, the derivatives of 
the M and C matrices with respect to the Q matrix. 


For the first step, it can be shown that 


dA aa. 
n 
where 
ee up, D S KINO SR D.D (d An + dAn) 

R : (ap, D + 5, 4D) 22D Ba- D DOEN B,- 

PL Drei (aD' Ba- +a B,1)) . 
and 

dB, = —~ LORD: DAD 
nr 

where 


dB ag = dB) S+B/ ds. 
S =D (up ` 


nyXxny 


m SN l 


dS = dÔ (up ` = SN +Ô up ` - Mm ` (ab) Do! dd), — hd, l. 


To compute dC and dM in the second step, we adapt the methodology de- 
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ul 
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scribed in Anderson, McGrattan, Hansen and Sargent (1996) and obtain 


vec(dC) = [(Q 81) +(1@CQy) + (CQ, 81) — (18 Qu)] 
— (I Q C) vec (dQ:1) — (C’ Gi C) vec (dQi2) 
+ vec (dQ21) + (C’ Q D vec (dQ22) 


and 
vec (dM) = vec (dQ11) + (C 8 T) vec (dQ12) + (I ® dQ12) vec (dC). 


A.5 Risk premia 
A.5.1 Holding premia 


We define the one-period holding premium en, on an n-period bond purchased at t 


as the expected holding return of that bond over one period, less the risk-free rate: 
ent = E; [In (PETE) — In (p?)] — re- 
Using the bond pricing equation, this can be written as 
Ge IS. rh = Ap, peens, e (LP, Du Z: 


A.5.2 Forward premia 


The one-period forward premium Yp at t for maturity n is defined as the differ- 
ence between the implied one-period forward rate n periods ahead, fn,t, less the 


corresponding expected one-period interest rate: 


Vnt = fnt — Er [rtn]. 
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The implied forward rate is given by 


Lu = In(p?)—In (pp) 


2 E = T 4 [Br DD (up ` = EA) $ A Ze 


while the expected short rate is 


The one-month forward premium is therefore 


Wnt 


fnt — Er [risa] 


os (NEES Sa 
(pc z T dh 


E — Dn (up ` S SU +n (7-DM"D"")| Z 


A.5.3 Yield risk premia 


The n-maturity yield premium at t, Aur, can be defined as the average of the forward 
premia up until t +n — 1, Le wnt = 4 KEREN ar, This is given by 


n—1 
n 2 
(B; - BÔ (MÔ - =) +4’ (7-DM'B")) Z]. (15) 
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Table 1: Parameter estimates 
(Sample period: Feb. 1975-Dec. 1998) 


Parameter Point estimate Standard error 


p 0.976 0.015 
B 2.087 0.855 
y 1.243 0.925 
Un 0.132 0.011 
dx X 10? 0.038 0.054 
Uz 0.303 0.029 
C; 0.027 0.023 
bat 0.999 — 
Ge X 10? 0.014 0.001 
ga X 10? 0.040 0.001 
Or X 10? 0.022 0.001 
Oz X 10? 0.097 0.004 
ot x 10? 0.019 0.011 
o% x 10? 0.025 0.014 
oi x 10? 0.023 0.001 
o? x 10? 0.023 0.001 
An —0.421 0.306 
o,2 —0.587 0.345 
do,3 4.431 2.565 
An —1.693 1.438 
Au x 1072 
ar? r T x 
art 0 0 0.976 0.912 
(-) (-) (0.264) (0.175) 
r | —35.354 18.955 -8.535 0 
(12.920) (7.184) (3.396) En 


m | 152.232 —88.432 41.470 0 
(56.674) (29.876) (12.488) (-) 


£ 0 2.094 0 2.559 
(-) (1.070) (-) (0.908) 


Standard errors in parentheses 


Asymptotic standard errors are based on the 
outer-product estimate of the information ma- 
trix. The estimates of the lag coefficients for 
inflation and output are not reported. 
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Table 2: Out-of-sample output and inflation forecasting performance: RMSEs 


Forecast horizon 


3-month 6-month 9-month 12-month 
variable HTV-M HTV HTV-M HTV HTV-M HTV HTV-M HTV 
x 1.242 1.508 1.776 2.416 2.064 2.889 2.154 2.980 
T 0.393 0.379 0.519 0.490 0.679 0.602 0.900 0.751 


RMSEs for out-of-sample forecasts between 1995:01 and 1998:12, based on parameter estimates for 
1975:02 - 1994:12. HTV-M denotes the macroeconomic model represented by equations (1)-(3) in the 
text (this model is estimated using the the Bundesbank’s price norm as the inflation target in the policy 
rule), and HTV denotes our structural macro model. 
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Table 3: Out-of-sample yield forecast performance: RMSEs 


1-month forecast horizon 


maturity RW VAR Ao (3) AP HTV 
1 month 0.148 0.182 0.151 0.146 0.129 
3 months 0.173 0.177 0.181 0.178 0.220 
1 year 0.194 0.211 0.319 0.271 0.270 
3 years 0.252 0.267 0.254 0.256 0.236 
7 years 0.220 0.237 0.331 0.320 0.384 


3-month forecast horizon 


maturity RW VAR Ao (3) AP HTV 
1 month 0.299 0.337 0.303 0.263 0.210 
3 months 0.345 0.358 0.426 0.446 0.292 
1 year 0.395 0.433 0.582 0.569 0.383 
3 years 0.448 0.544 0.462 0.475 0.397 
7 years 0.379 0.454 0.428 0.465 0.447 


6-month forecast horizon 


maturity RW VAR Ao (3) AP HTV 
1 month 0.458 0.425 0.452 0.410 0.296 
3 months 0.512 0.477 0.652 0.743 0.404. 
1 year 0.574 0.604 0.829 0.873 0.529 
3 years 0.624 0.765 0.669 0.885 0.534 
7 years 0.521 0.684 0.521 0.706 0.479 


9-month forecast horizon 


maturity RW VAR Ag (3) AP HTV 
1 month 0.610 0.548 0.588 0.523 0.432 
3 months 0.666 0.656 0.828 0.975 0.559 
1 year 0.733 0.786 1.012 1.138 0.689 
3 years 0.782 0.986 0.873 0.950 0.678 
7 years 0.678 0.916 0.719 0.848 0.580 


12-month forecast horizon 


maturity RW VAR Ao (3) AP HTV 
1 month 0.747 0.678 0.706 0.608 0.604 
3 months 0.793 0.857 0.956 1.154 0.730 
1 year 0.854 1.002 1.140 1.333 0.829 
3 years 0.842 1.198 1.002 1.139 0.744 
7 years 0.806 1.184 0.920 1.091 0.680 


RMSEs for out-of-sample forecasts between 1995:01 and 1998:12, based on parameter estimates for 
1975:02 - 1994:12. "RW" are random walk forecasts, "VAR" is an unrestricted VAR(3) including 
the same variables as our model, "Ao (3)" is a canonical essentially affine Gaussian three-factor model, 
"AP" denotes the Ang-Piazzesi (2003) Macro Model (estimated using our macro data, but with inflation 
expressed in y-o-y terms), and "HTV" denotes our structural macro model. 


ECB 
2 Working Paper Series No. 405 
November 2004 


Table 4: Out-of-sample yield forecast performance: Trace MSEs 


forecast horizon 
(months) 


SEAN anrkRwne 


11 
12 


RW 


(E 
16.4 
25.5 
33.6 
41.0 
52.7 
64.0 
74.4 
87.3 
99.9 

109.6 
117.9 


VAR 


8.5 
19.2 
33.5 
44.4 
55.7 
65.7 
81.3 
94.2 

113.8 
135.8 
159.4 
181.2 


Ao (3) 


11.9 
23.8 
36.3 
47.4 
58.1 
73.2 
88.3 
102.8 
120.0 
136.7 
150.9 
164.2 


AP 


10.6 
23.1 
37.2 
51.1 
65.7 
85.6 
105.7 
125.9 
149.1 
171.9 
193.0 
214.9 


HTV 


12.3 
17.7 
22.8 
26.3 
29.8 
37.6 
45.5 
52.9 
63.7 
75.0 
84.4 
93.6 


The trace MSE statistics are for out-of-sample forecasts between 1995:01 and 1998:12, based on pa- 
rameter estimates for 1975:02 - 1994:12. "RW" are random walk forecasts, "VAR" is an unrestricted 
VAR(3) including the same variables as our model, "Apo (3)" is a canonical essentially affine Gaussian 
three-factor model, "AP" denotes the Ang-Piazzesi (2003) Macro Model (estimated using our macro 


data, but with inflation expressed in y-o-y terms), and "HTV" denotes our structural macro model. 
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Table 5: Tests for superior out-of-sample predictive ability of yield 
forecasts from the HTV model compared to four different bench- 
marks 
1 month forecast horizon 
maturity RW VAR Ao (3) AP 
1 month 0.039 0.001 0.036 0.104 
3 months 0.995 0.995 0.942 0.963 
1 year 0.999 0.999 0.007 0.207 
3 years 0.002 0.053 0.002 0.003 
7 years 0.997 0.991 0.905 0.790 
3 month forecast horizon 
maturity RW VAR Ao (3) AP 
1 month 0.021 0.002 0.018 0.095 
3 months 0.022 0.135 0.003 0.002 
1 year 0.237 0.397 0.000 0.000 
3 years 0.002 0.001 0.002 0.002 
7 years 0.799 0.554 0.663 0.461 
6 month forecast horizon 
maturity RW VAR Ao (3) AP 
1 month 0.035 0.102 0.036 0.078 
3 months 0.017 0.322 0.005 0.002 
1 year 0.004 0.379 0.001 0.001 
3 years 0.003 0.001 0.002 0.001 
7 years 0.133 0.044 0.246 0.189 
9 month forecast horizon 
maturity RW VAR Ao (3) AP 
1 month 0.038 0.201 0.043 0.105 
3 months 0.024 0.296 0.006 0.001 
1 year 0.015 0.215 0.002 0.000 
3 years 0.003 0.000 0.000 0.000 
7 years 0.024 0.001 0.024 0.044 
12 month forecast horizon 
maturity RW VAR Ao (3) AP 
1 month 0.049 0.323 0.074 0.391 
3 months 0.086 0.250 0.012 0.001 
1 year 0.115 0.112 0.003 0.000 
3 years 0.006 0.001 0.000 0.000 
7 years 0.009 0.000 0.004 0.010 
The table shows p-values corresponding to test statistics for superior forecast ability of the HTV model, 
compared to each of the four different benchmarks listed in the tables, calculated according to White’s 
(2000) "reality check." We use a squared forecast error loss function when implementing the test. The 
null hypothesis is that the expected differential between the forecast loss of the benchmark and that of 
the HTV model is smaller than or equal to zero. The p-values are based on the stationary bootstrap 
approach, with 50,000 resamples of the loss differential series (using a smoothing parameter of 1/12). 
Bold figures denote rejection of the null at the 5 percent level. 
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Figure 1: Data used in the estimations 
(a) Macro data 


15 —— 7 7 


Inflation 


— — Output gap 


—10 


1979 


1987 


1995 


-15 


1975 1983 1991 1999 
The inflation and output gap series have been multiplied by 100. 


The sample period is January 1975 to December 1998. 
(b) Yield data 
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German term structure data over the sample period January 1975 
to December 1998 (percent per year). 
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Figure 2: Estimated inflation target and announced Bundesbank price norm 
6 r r 7 
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nounced upper and lower bounds for the price norm, an average 
of these is shown in the figure. 


Figure 3: Impulse responses from inflation target shock 


Output gap nflation 


All responses are expressed in percentage terms. The inflation 
and short rate responses are expressed in annual terms. Dotted 
lines are 95% confidence bands. The inflation target was shocked 
by one standard deviation (around 0.2% p.a.). 
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Figure 4: Impulse responses from monetary policy shock 


All responses are expressed in percentage terms. The inflation 
and short rate responses are expressed in annual terms. Dotted 
lines are 95% confidence bands. The short-term interest rate was 
shocked by one standard deviation (around 48 bps.). 
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Figure 5: Impulse responses from inflation shock 


Oulpul gap nflation 


All responses are expressed in percentage terms. The inflation and 
short rate responses are expressed in annual terms. Dotted lines 
are 95% confidence bands. Inflation was shocked by one standard 
deviation (around 0.26% p.a.). 
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Figure 6: Impulse responses from output shock 


Oulpul gap nfla 


All responses are expressed in percentage terms. The inflation 
and short rate responses are expressed in annual terms. Dotted 
lines are 95% confidence bands. The output gap was shocked by 
one standard deviation (around 1.2%). 
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Figure 7: Initial response of yield premia to macro shocks 
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The figure shows the one-month ahead response of the yield pre- 
mia Gin, at maturities n up to 84 months, to one standard devia- 
tion shocks to the four macro factors. The premia are expressed 


in annual percentage terms. 
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Figure 8: Estimated yield premia and components of premia 


De—meaned 1-year premium De—meaned 1-year premium 


interest rate and target 


P] Hei 

N i 1 1 1 N 1 1 1 

| 1975 1979 1983 28 1991 1995 1 11975 1979 1983 1987 1991 1995 999 
De- meaned year premium De-meaned 7-year premium 


1991 1995 


1983 1987 1991 1995 1 ) 11975 1979 1983 


| 1975 1£ 979 98¢ 
The solid lines are the estimated (de-meaned) yield premiums w,, during the sample period, 
for maturities n = 12 and 84 months, expressed in annual percentage terms. The dashed 
lines show the portions of the premia that are due to selected macro factors or combinations 
of such factors. 
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Figure 9: Model-implied CS projection coefficients: ” LPY (i)” 
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Empirical estimates of the CS long-rate coefficients d. in ZEN — 
yy = Pn (Y? — rt) / (n — 1), plus corresponding model-implied cof- 
ficient values. The "population" coefficients are the theoretical 
values based on our estimates; the MC coefficients are the mean 
estimates from 1000 series of the same size as the sample, simu- 
lated from our model. The bands around the MC mean estimates 
are 5% confidence bands. 
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Figure 10: Model-implied risk-premium adjusted CS coefficients: ” LPY (ii)” 
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The figure shows the estimates of the Campbell and Shiller 

(1991) long-rate coefficients ¢, in the regression ey -y = 

dn (yf — re) / (n — 1) for our sample, along with the correspond- 

ing risk-premium adjusted model-implied cofficient values based 

on our parameter estimates. 
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